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Abstract. We derive and solve a Boltzmann equation governing the cosmological evolu- 
tion of the number density of current carrying cosmic string loops, whose centrifugally 
supported equilibrium configurations are also referred to as vortons. The phase space is 
three-dimensional and consists of the time variable, the loop size, and a conserved quantum 
number. Our approach includes gravitational wave emission, a possibly finite lifetime for the 
vortons and works with any initial loop distribution and for any loop production function. 
We then show how our results generalize previous approaches on the vorton excess problem 
by tracking down the time evolution of the various sub-populations of current-carrying loops 
in a string network. 
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1 Introduction 

Cosmic strings are one of the expected, and actively searched, signatures of early universe 
physics in cosmology [1-12]. Originally, they were considered as line-like topological defects 
produced during the phase transitions associated with the spontaneous breakdown of Grand 
Unified symmetries [13, 14]. More recently, it was realized that very similar objects could 
also be generated within Superstrings Theory and formed at the end of brane inflation [15— 
23]. In their simplest version, cosmic strings exhibit Lorentz invariance along the string 
worldsheet and the two prototypical examples are the Abelian-Higgs string and the Nambu- 
Goto string. Both of these models are simple enough to allow for numerical simulations of 
cosmic strings evolution in Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetimes [24- 
32]. These simulations have shown that, due to their so-called scaling regime, long strings 
never dominate the energy density of the Universe and evolve as oc 1/i 2 , instead of 
the naively expected 1/a 2 behaviour for one-dimensional objects (t being the cosmic time 
and a the scale factor). In this regime, the unique model parameter is the string energy 
density per unit length, U, which is constrained by Cosmic Microwave Background (CMB) 
measurements. For the Abelian-Higgs string network, Ref. [33] reports a two-sigma upper 
bound G N U < 4.2 x 10~ 7 , C N being the Newton constant. For Nambu-Goto strings, the 
scaling evolution takes place via the emission of loops whereas the Abelian-Higgs string 
scaling seems to mainly rely on particle emission [30, 32]. In addition to be a priori model- 
dependent, the expected cosmological distribution of loops is a matter of debate. As shown 
in Refs. [31, 34], a copious number of small loops seen in Nambu-Goto simulations can be 
linked to some relaxation effects coming from the initial conditions, and as such they are 
not expected to persist over cosmological timescales. However, in addition to these transient 
loops, Ref. [31] was the first to exhibit a sub-population of loops in scaling, i.e having an 
energy density evolution similar to the one of the long strings. Their distribution was also 
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found to match the analytic and independent findings of Refs. [35-37]. Contrary to the 
former, these loops are expected to persist during the cosmological expansion. Those results 
were recovered in Refs. [38-40], albeit with some differences in the distribution, while the 
authors of a more recent simulation presented in Ref. [41] arguably claim that all loops but 
the Hubble-sized ones (the so-called Kibble loops) should ultimately disappear 1 (see Ref. [9] 
and references therein for a review). 

Even if we could infer the exact loop production function from the string simulations, 
free of any numerical uncertainties, one would have to keep in mind that the dynamics of 
the smallest length scales still involves physical processes that are not included in numerical 
simulations. For instance, gravitational wave emission and gravitational backreaction are 
crucial ingredients which require the use of analytic and semi-analytic models to be dealt 
with [42-49]. In particular, based on the loop production function of Refs. [35, 36], the 
Boltzmann approach of Ref. [37] has been extended in Ref. [46] to include the effects of the 
initial conditions and gravitational backreaction into the evolution of Nambu-Goto cosmic 
string loops. A Boltzmann treatment for loops consists in describing their evolution in phase 
space, e.g. (t, £) if I denotes the loop size, in order to extract their number density distribution 
n(t,£). Such a method has first been used in Ref. [50] to study the generation of long strings 
from networks having only loops. 

In this paper, we go one step further and present a new Boltzmann treatment dealing 
with the so-called current carrying cosmic string loops, a class of uttermost importance for 
cosmology [51]. These objects naturally appear as soon as the strings exhibit coupling with 
bosonic or fermionic particles [52-57] . This results into the breaking of Lorentz invariance 
along the worldsheet thereby allowing the existence of centrifugally supported loops called 
vortons [58, 59]. If these objects are stable over cosmological timescales [60-63] they could 
come to dominate the energy budget of the Universe and this has been used in Refs. [64, 65] 
to severely constrain the energy scale at which they may be formed. In these papers, it was 
nevertheless assumed that most of the loops, which can potentially become vortons at later 
times, are formed once and for all at the string forming time. In view of the new numerical 
results described above, one may wonder whether more drastic constraints are not in place 
when one takes into account the loops produced during the cosmological evolution of the 
string network. 

In the following, we tackle this problem and solve a Boltzmann equation describing the 
production and decay of current carrying loops. Our approach includes gravitational radia- 
tion, finite life-time vortons with any initial distributions and any loop production functions. 
After having presented our working hypothesis and written down the Boltzmann equation, 
we first present the Liouville solution in Sec. 2.3, i.e. without a loop production function. 
This allows us to physically discuss the different population of loops present at any time. 
In Sec. 2.4, we solve the complete Boltzmann equation, in presence of a loop production 
function. This is our main result and it is presented in Eqs. (2.24), (2.37) and (2.39). In 
Sec. 3, we generalize the results of Ref. [64] to the case of perfectly stable vortons created 
from any initial loop distribution. This allows us to track down the time evolution of the 
so-called proto- vortons, doomed loops and vortons. We conclude in the last section. 

1 The string simulation of Ref. [41] has currently the longest dynamical time, but at the expense of having 
the poorest resolution. Their discretized strings having only two points per correlation length as opposed to 
twenty for Ref. [31] and hundreds for Ref. [39]. Among others, a low resolution has precisely the effect of 
artificially cutting the loop distribution at small scales (see Fig. 2 in Ref. [31]). 
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2 Boltzmann approach 



2.1 Assumptions 

Our main assumption is that the loop production mechanism remains similar to the Nambu- 
Goto case even though the strings are now current-carrying. This can be motivated from 
the macroscopic covariant formalism of Carter [52, 66] as the existence of a current can be 
viewed as modifying the Nambu-Goto relation U = T (where T is the string tension) into 
a more general equation of state U(T). Provided the currents are not too strong, one has 
U ~ T [53, 67] and the string dynamics remain similar to the Nambu-Goto one. Let us 
notice that some cosmic superstring models clearly violate this assumption, as for instance 
those developing Y-j unctions [47, 68-70]. 

Current-carrying loops are assumed to be created with two conserved quantum and 
topological numbers N and Z associated with the existence of a conserved current flowing 
along the string [71] . Following Ref . [64] , we assume Z ~ N with 



Here £ is the physical length of a loop at formation and A the typical wavelength of the 
carrier field fluctuations, which is considered to be constant (originating from microphysics, 
it can be assumed to be of the order of the Compton wavelength of the condensed particle). 
This expression should be valid provided A remains much smaller than the mean interstring 
distance. In this paper, we will make no assumptions on the explicit shape of the loop 
production function and simply call it V(£,t). As a well-motivated simple example, one can 
pick that of Refs. [35-37] for small loops, adding a Dirac distribution centered at the mean 
inter-string distance for the Kibble loops. 

2.2 Evolution in phase space 

Denoting by n(£, N, t) the number density distribution of cosmic string loops of size £, and 
conserved quantum number ./V at cosmic time t, we can write 



where V(£,t) is the loop production function discussed above, which is a collision term from 
the Boltzmann point of view, a? accounts for space-time expansion while J~(£,t) encodes 
any additional phase space distortion. The Dirac distribution ensures that all loops freshly 
formed carry the expected quantum numbers provided by Eq. (2.1). Due to gravitational 
radiation, a current-carrying loop of size £{t) shrinks at a constant rate until it eventually 
becomes a vorton, i.e. a state centrifugally supported by its current. Its fate depends on the 
microscopic model under scrutiny. For the sake of generality, we are assuming that vortons 
can decay by an unspecified mechanism at a constant rate 7 V , the completely stable situation 
corresponding to 7 V — > 0. As a result, we model the evolution of £{t) for a loop of conserved 
number N by 



where ®{x) denotes the Heaviside function. The gravitational wave emission rate, 7d, is 
assumed to be the same as for Nambu-Goto strings 






(2.2) 




(2.3) 



7d ^ TG N U, 



(2.4) 
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where T is a constant of order 10 2 which depends on the string structure, geometry and 
dynamics [72, 73]. Having in mind the limit 7 v /7d — * for the actual vorton situation, one 
can safely assume 7 v /7d <C 1 in numerical evaluations. The quantity £ V (N) is the physical 
length of vortons typically given by [64, 65] 



tv(N) 



NZ 



N 



U 



(2.5) 



In the phase space (t,£,N), Eq. (2.3) does not preserve volume. All loops contained in a 
given interval A£{ n [ evolve according to Eq. (2.3) into a smaller A£(t) interval once £(t) < £ v . 
We are therefore in presence of a generalized Liouville evolution and J'(£,t) is the Jacobian 
of the transformation mapping £- mi to £(t) = £,(£- m i,t), i.e. 



J(£,t) 



(2.6) 



The function £(£im,t) is the characteristic curve starting at £- 1Ti i and solution of Eq. (2.3) (see 
Fig. 1). Concerning the conserved number N, by definition we have dN(t)/dt = along 
any loop trajectory and there is no phase space volume distortion along the N direction. 
Combining Eqs. (2.2) and (2.3) gives the Boltzmann equation 



0_ 
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a 3 J(£,t) 
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d£dN 
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+ 7v6 



a 3 J(£,t)V(£,t)5[ N 
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a 3 J(£,t) 
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d£dN 



(2.7) 



In order to make contact with previous works for Nambu-Goto strings [46], one can also 
express the loop sizes £ in units of the cosmic time t. For this purpose, we will also use in 
the following the variable 7 and the density function T defined as 



7 (e,t) 



T = J 



d 2 n 
d£dN 



(21 



For a given value of N, Eq. (2.7) exhibits two typical length scales, one is £ V (N) already 
given by Eq. (2.5) and the other is 



£ p (N) = XN 2 



(2.9) 



which gives the size of a loop with conserved number N at production time. For N = 0, we 
have £ v = £ p = 0, otherwise 

U(N) « £ p (N) , (2.10) 

assuming y/U 3> A -1 . In the following, we first discuss the Liouville equation, i.e. Eq. (2.7) 
without the production term, before solving the full Boltzmann equation. 

2.3 Liouville evolution 

In order to get intuition on the solution, we solve first the Liouville equation, i.e. without 
including the loop production. Under this assumption, one should recover the results of 
Ref. [64] since all the produced vortons come only from the initial loop distribution (see also 
Sec. 3). 
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'-ini £ 

Figure 1. Characteristics £(i?i n i, t) of the Liouville operator at fixed N in the plane (i, I). At any time 
t, the Cauchy problem is solved by specifying initial conditions on the surface t = and propagating 
them along the characteristics in the "+" and "— " regions. Notice that the solution in the "0" region 
is determined by specifying initial condition on the line £ = £ v , which can still be deduced from t = t ml 
after evolution through the "+" domain. 



2.3.1 Jacobian 

In order to determine J(t,t), the function £(4m,*) is needed. Integrating Eq. (2.3) gives 



7v*) + ®(4n " U) { (4n " 7d*) ©(Ani ~ ly - 7d*) 



7v 



4 " 7vi + — (Ani - ty) 

7d 



e( 7d t + 4 



(2.11) 



from which we deduce that 



{£- mi ,t) = e(£ v -£ ini ) + @{£ illi -£ v ) 



G(4 



7d*) + —Q(idt + 
7d 



(2.12) 

where we have set ti n ; = 0. 

The Jacobian J{£, t) is obtained by taking the modulus of the above equation, up 



to a change of variables (£i n ut) 



?, t). This requires the inversion of Eq. (2.11) to get 



4 



£ 1 (£, i). After some algebra, the final expression simplifies considerably and reads 



J(£,t) = Q(£-£ v ) + Q(£ v 



6(4 - I - 7v i) + — &(£ + l v t - £ v ) 
7d 



(2.13) 



which merely reflects the changes of slopes of the decaying processes (see Fig. 1 discussed 
below). At any time t and for any value of N, the Jacobian is unity for all loop sizes satisfying 
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£ > £ y (N). It then switches to the value 7 v /7d in the domain £ v — 7 v i < £ < £ v and returns to 
one again for £ < £ v — j v t. Fig. 1 shows these three domains in the (t, £) plane, respectively 
denoted as "+", "0" and Their origin lies in the various histories a £; n i-sized loop 

starting at t = t- in i may have. From Eq. (2.11), all loops starting with £- m i < £ v ("— " domain) 
decay with a rate given by 7 V as being already centrifugally supported. The loop number 
density distribution per A£ interval is thus preserved and the Jacobian is unity. At any time 
t, all loops having £ > £ v are necessarily coming from £- m \ > £ v . Since they never crossed £ v , 
they have shrunk only by gravitational radiation, again at a constant rate given by 7d, the 
Jacobian is also equal to unity. Only the loops in the "0" domain have a two steps history. 
They started with £- m i > £ v , then radiated gravitational waves until they became vortons 
by hitting £ = £ v . Once vortonized, their energy emission rate being j v <C 7d, they start 
accumulating in this region and the Jacobian takes the value 7 v /7d ^ 1- 

2.3.2 Loop number density distribution 

We are now in the position to solve Eq. (2.7) without the collision term, i.e. in absence of 
loop production. This is a first order partial differential equation whose characteristics are 
given by Eq. (2.11). 

According to the values of £, it is convenient to introduce the new coordinates (r, x) 
and (r, y) defined as 

V£>£ V (N) => r = t, x = ld t + £, 
V£<£ V (N) =► T = t, y = j v t + £. 

In terms of these variables, the Liouville equation, i.e. Eq. (2.7) with V(£,t) = 0, takes the 
simple form 

8 <^U . (2.15) 



The solutions are 



dr 

a\t)F{t,£> £ V ,N) = I(£ + 7d t,N), 



( 2 - 16 ) 



a 6 (t)F{t, £ < £ v , N) = J(£ + 7v t, N), 



where /(•,•) and «/(•,•) are two arbitrary functions still to be determined. Given the loop 
distribution at some initial time £i n i, 



2 



71 

(t ini ,£,N)=Mm(£,N), (2.17) 



d£dN 

one can determine these two functions in their respective domain of definition. Plugging 
Eq. (2.16) into Eq. (2.17) and making use of Eqs. (2.8) and (2.13), all evaluated at t = t- m \ 
gives 

I(z > z ini , N) = af ni AA ini (2 - 7diini, N), (2.18) 

where z- m \ = £ v +7diini- This expression completely fixes the solution inside the "+" domain 

as 

a 3 (t)g^(t,£,N) = al i M ini [£ + ld (t-t ini ),N], 

£ > £ V (N). 

According to the discussion of the previous section, the solution in the "0" domain is uniquely 
determined by propagating this solution to £ = £ y , at all times. The matching condition 
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Figure 2. Liouville time evolution of the loop distribution given by Eqs. (2.19), (2.22) and (2.24) with 
7d = 1. At t = Uui, we assume a gaussian initial distribution centered at £ = £ p = 5 with a standard 
deviation of a p = 1 (full black line). We also assumed perfectly stable vortons, so that j v = 0, and 
we smoothed the step function appearing in Eq.(2.3) with Q(£ — £ v ) — > {1 + tanh [(£ — £ V )/X]} /2, 
and A = 0.5: this is the reason why the distribution is non-vanishing for values of £ smaller than 
£ y . As time evolves, the distribution is shifted to the left, where it starts rapidly feeling the vorton 
characteristic scale £ v , here numerically set to 2. Each curve represents the distribution one unit 
of time after the previous curve, i.e., in the units used, the initial gaussian is for t = 0, and then 
t = 1, 2, . . . , 8 (orange dotted, pink dashed, and so on, respectively). As time passes, the distribution 
evolves to a very sharp peak located around £ v , which gives a Dirac distribution 5(£ ~ £ v ) in the limit 
A -> 0. 



implies 2 
such that 



afniMni^v + 7d(* - tird), N] = J(£ v + 7v t, N), 



J(z > z v ,N) = af ni Mm 



7v 



(2.20) 



(2.21) 



Here, we have defined z N = ^ v +7viini- Plugging this expression back into Eq. (2.16), restricted 



The conserved quantity is a T which already includes the Jacobian factor, see Eq. (2.15). 
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on I + -y v t > z v (and £ < £ v ) completely fixes the solution inside the "0" domain 



didN^ ' ' 1111 7v 



,y + — (£-£y)+J d (t-t illi ),N 

7v J (2.22) 

£ V (N) - lv (t - t ini ) < £ < £ V (N) . 



"-" domain for which the initial conditions are fixed at t = t- m i and 



At that point, J(z,N) is not yet determined for z < z v . However, this part is exactly the 

ai 

J{z < zy, N) = a? ni N ini (z - 7v t ini , N). (2.23) 

The solution finally reads 

a3{t) dMN it,e,N) = a ?nrMni[^ + 7v(t " t ini ),N] , 

£ < £ V (N) - jy(t - tini). 

The loop distribution at all times t > t; n i, length £ and conserved quantum number N is 
then given by the set of Eqs. (2.19), (2.22) and (2.24). Although the solutions in region "+" 
and "-" could have been intuitively guessed as they only shift the initial loop distribution 
according to the gravitational and vortonic energy emission rates, Eq. (2.22) shows that loops 
tend to accumulate into the region "0" . 

This is exemplified on Fig. 2, showing the time evolution of an initially gaussian distri- 
bution and producing perfectly stable vortons, case to which we now turn. 

2.3.3 Perfectly stable vortons 

In the limit of perfectly stable vortons, i.e. 7v — > 0, Eq. (2.22) is singular and the loop density 
distribution becomes infinite but into a phase space region which becomes infinitely small 
(see Fig. 1). This is not problematic as the physical quantity of interest is the total number 
of loop in that region. We can therefore define the number of "vortons" having a charge N, 
at time t by 

f iv d 2 n 

vut ' N) -L,,-^ {tXN)M ' (2 - 25) 

where "ini" is a reminder that, in absence of loop production, they can only come from the 
initial loop distribution. Using Eq. (2.22), defining for each time slice the new variable 

a = ^(£-£ v )+j d (t-t iQi ), (2.26) 
7v 



the integral reduces to 

,3 /•Td(t-W) 

a 3 (t) 



V ini (i, AO = -^sL / A/- ini (4 + a, N)da , (2.27) 



which is well-defined and finite. For an homogeneous initial distribution (A/i n i = C tc ) one 
recovers the intuitive result that the number of vortons grows linearly with time as a 3 Vmi oc 
7d(* - £ini)- 



- 8 - 



We conclude that, in the limit of perfectly stable vortons 7 V — > 0, the loop density 
distribution reads, for all I, 

hm a 3 (t)^(t,£,N) = Q(£ v - £) af ni Mni(^, N) 

+ @(e - e v ) af ni M ini [e + -y d (t- t ini ) , n] (2.28) 

/•7d(*-*ini) 

+ S(£ - £ v ) af ni / Mni(4 + a, N)da . 



Jo 

2.4 Boltzmann evolution with loop production 

We now solve Eq. (2.7) in full generality by separating the phase space into the same three 
domains. 

2.4.1 Formal solution 

For all £ > £ V (N), Eq. (2.7) can be recast in terms of the variable (r, x) introduced in 
Eq. (2.14) as 



a A (r)V{x - 7d r, r)5 1 N - ^ — ^_ I , (2.29) 

the Jacobian being unity there. The delta function has an argument depending on r such 
that it simplifies to 



5 [N-r-^ =%^V-r p (x,iV)], (2.30) 



where £ P {N) is given by Eq. (2.9) and 



x-£JN) 

t p (x,N) = — — - . 2.31 

7d 



The loop distribution in the "+" domain is then solution of 



2a 3 [r p (x, N)] V[£ p (N),r p (x, N)} 8[r - r p (x, N)] , (2.32) 

and reads 

alF(r, x, N) = I(x, N) + 2a 3 [T p (x, N)] V[£ P (N), r p (x, N)] 6[r - t p (x, N)} . (2.33) 

As for the Liouville equation, I(x, N) is a unknown function that has to be determined from 
the initial conditions. In terms of the original variables, the previous equation reads 

:(t,£,N)=I(£ + ld t,N) 



d£dN 



+ it a .( t+ £z4W < ., t+ £z4W-<), 

7d^ V 7d / V 7d / 

£ > £ V (N). 
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For all £ < £ V (N), we can similarly use the new variables (r, y) denned in Eq. (2.14) to 
solve Eq. (2.7). In fact, since £ p 3> £ v (or null), the source term always vanish in that domain 
and the solution is the same as in Sec. 2.3, i.e. 



a 3 J(£,t) 



ff 2 



n 



dtdN 



(t,£,N) = J(£ + lv t,N), 
£ < £ V (N), 



(2.35) 



where the function J(z, N) has again to be determined from the initial conditions. Notice 
that the Jacobian J{£,€) is either unity or equal to 7 v /7d according to the domain "-" or 
"0", respectively. The formal solution is therefore given by both Eqs. (2.34) and (2.35). 



2.4.2 Full solution from specified initial conditions 

As for the Liouville equation, we assume the loop distribution to be known at some initial 
time tim and given by Eq. (2.17). 

For all £ > £ v , i.e. in the "+" domain, Eq. (2.34) evaluated at t = t; n i fixes the function 



I(z > z ini ,N) = af ni A/"ini(z-7dtini, A0- 



21 



p a 3 



z_-£ 1 
7d 



x9 t ini + 



£ P -z 

7d 



V 



z-L 



7d 



(2.36) 

where z- iri ^ = £ v + 7dimi- Notice that the Heaviside function ensures that the argument of the 
scale factor a(-) is always larger or equal than ti n ;. Plugging I(z,N), back into Eq. (2.34) 
gives the wanted result 



a\t) 



2 



d£dN 



n 2£ 

(t, £, N) = af ni Mm[£ + 7d(* - t^N] + 1 t 



7d 



xV[£ p) t 



7d 



{e(£ p -£)- Q[£ p - £ - 7d (t - t ini )] }, 

£ > £ V (N). 



(2.37) 



For all £ < £ V (N), we start from the solution (2.35). However, as in Sec. 2.3, there are 
two sub-cases corresponding to the domain "0" and "-" of Fig. 1. 

On one hand, for £ v — j v (t — ti a i) < £ < £ v , we determine the J function by matching 
Eq. (2.35) to Eq. (2.37) at t = £ v and for all times. One gets 



J(z > z v ,N) = af ni AA ir 



+ ^( z -£ v )- ld t inh N 

7v 



+ 



2f 



p a 3 



7diV 



xV £ 



z-L 



-P' 



7v 



£p £\ 
id 







£p £\ 
id 



7d 



z-£ v 
7v 

z_-£ 1 

7v 



£ — £ 
7d 



(2.38) 



where, as before, z v = ^v + 7v^ini- This completely fixes the solution in the "0" domain which 
reads 



7v^ 



7d ; 



£,)+ldit-t inl ) ) N 



7v 



7d 



V 



■pi 



7v 



7d 



(2.39) 



x 0(£ P - £ v ) - 9 



+ ^(£ v -£)- ld (t-t ini ) 

7v 

£ V (N) - 7v (t - t ini ) < £ < £ V (N). 
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Notice that the Heaviside function Q(£ p — £ v ) = 1 under our hypothesis t p 3> t v - 

On the other hand, when t < £ v — 7 v (i — i; n i), the function J(z < z v , N) is set by the 
initial loop distribution at t = t{ n [. Since the source terms vanish in that region, combining 
Eqs. (2.17) and (2.35) yields exactly Eq. (2.23) for J. As a result, the loop distribution in 
that domain is the same as the one derived from the Liouville operator in Eq. (2.24). 

To summarize, starting from the initial distribution of Eq. (2.17), in presence of a loop 
production function, the loop number density distribution at time t, for all £ and N is given 
by Eqs. (2.37), (2.39) and (2.24). 

2.4.3 Perfectly stable vortons with loops production 

As for the Liouville solution, our results also apply to the particular limit 7 V — > 0. The 
number of vortons can be defined as in Eq. (2.25) and is obtained by integrating Eq. (2.39) 
over all loop sizes intercepting the "0" domain at time t. The first term is the same as in 
Eq. (2.27) while there is an additional contribution coming from the loop production function. 
Defining 

t p (N)-U.N) iV / 1 



the time required for a loop of size £ P (N) to reach £ V (N) by gravitational radiation, and the 
new variable 

t — f 

, = L^L +t , (2.41) 
7v 



one gets 

21 



V(t,N) = V ini (t,N) 



p 



a 3 (t)N 



! a 3 {^-t pv )r(£ p ,q-t pv )@^-t pv )d<i. (2.42) 



This expression does no longer depend on j v which makes it valid even for j v — > 0. It shows 
that, in addition to the vortons coming from the initial loop distribution, loop production 
incessantly feeds the vortons reservoir with a time delay given by t pv (N). 

The loop distribution when 7 V — > 0, in presence of a loop production function, finally 
reads for all £: 



7v->o o£oJS 



mr 
,3 



+ 0(1 - £ v ) af ni Mni[^ + 7d (t - iini) , N] 
+ @(£-£ v )^-a 3 (t- e -^-)v(£ p ,t 



7dN V 7d J V 7d 
x {e(£ p -£)- Q[l p - I - 7d (t - t ini )] } (2 ' 43) 

/•7d(t-*ini) 

+ 5(£- £ v ) af ni / M ini (£ v + <r, N)da 







2£ /"* 

+ / a 3 ( ? -t pv )P(£ p ,?-t pv )G( ? -t pv )d ? . 



Note that the initial loop distribution is already responsible for the so-called vorton excess 
problem, and we are adding a possibly larger term to it: as expected, looking in details at the 
loop distribution evolution constrains the models even more. We now turn to this question 
in the relevant cosmological framework. 
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3 Application to the vorton excess problem 



In Ref. [64], the relic abundance of vortons was estimated by assuming that loops are first 
formed during the string forming phase transition whereas they develop a current at a later 
time. Under some reasonable assumptions, Ref. [64] approximates the number density of 
vortons by n v oc nf(df/a) 3 with ri{ oc denoting the vorton abundance at formation 

if their distribution peaks around a particular length scale Lf. As we show below, our 
Boltzmann treatment allows to derive the complete current-carrying loop distribution without 
making assumption on their initial distribution. 

Under the same assumptions as in Ref. [64], we assume that loops are formed at an 
energy scale T x ~ y/U whereas the current-carrier condensation occurs at the lower temper- 
ature T cur ~ A" 1 . From the results of Sec. 2, we can readily write down the loop distribution 
at the time the strings become superconducting, i.e. at t = t CUT : 

AT cur (£) = -^L_M ini [£ + 7d (i cur - t iBi )], (3.1) 

where 

ATa = ^(l,t a ), (3.2) 

is the loop density distribution at t = t a . In Eq. (3.1), 7 d describes the loop gravitational 
decay, complemented or not by any other evaporation effects the loops may experience in the 
friction dominated regime. For simplicity, we will be keeping an unique 7 d in the following. At 
t = t CUT , current-carrier condensation occurs and those loops are endowed with the conserved 
number N induced by thermal fluctuations of wavelength A = 1/T cur : 



M cni (£,N)=M cm (£)5 [N-J- . (3.3) 




For £ > A, the loops can potentially become stable vortons. One should nevertheless require 
that £ V (N) = N/VU > A, i.e. 

N>N* = \VU. (3.4) 

In the opposite situation, N < N*, one may expect the current to disappear by some quantum 
processes and these loops are referred to as "doomed" . This effect can be phenomenologically 
included within our interpretation by assigning a vanishing effective vorton length for those 
loops, i.e. by imposing £ V (N < N*) = 0. For the large enough loops, following the hypothesis 
of Ref. [64], there is no loop production and the vortons having £ v > A are assumed to be 
perfectly stable 7 V — > 0. Under these assumptions, Eq. (2.43) becomes 



o2 



N 



+ 7d(£ - tcur) 

A 



/■7d(*-*cur) _ 

+ af ni <5[^ - £ V (N)} / A7 ini [£ V (N) +a + 7d (i cur - i ini )] 6 



n 



da. 



(3.5) 

Integrating this expression over N gives the number density of proto-vortons plus doomed 
loops (first line) and vortons (second line) of size £ at time t: 



• dn, „ f°° , d 2 
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Let us start by the number density of vortons. The first Dirac function in the second line of 
Eq. (3.5) allows to explicitly perform the integration over N. One has 



s[t - £ V (N)] = Vus(n- eVu^j 



(3.7) 



provided N > A*. In the opposite situation, £ V (N) = and the integral vanishes. All in all, 
one can use the above expression provided £\f\J > A*, i.e. for £ > A as expected. The vorton 
term now reads 



r1d(t— tour) _ / // i a 

af ni e(£ -X)VU J A7 ini [£ + a + 7d (t - t cur )} 5 ( ty/U - ' 



A 



2af ni ^ 2 ^ Q(t ~ A) e[£ T (t) - £} Mr 



A 2 y +7d(i-i CU r) 



(3.8) 



where £ T (t) is an increasing function of time defined by 



A 



2N? 



l+\l + 4N* 



7d{t - t cm ) 

A 



(3.9) 



As can be checked in Eq. (3.8), the remaining Dirac function yields a non- vanishing integral 
provided <7o = £{XU£ — 1) belongs to the interval [0,7d(t — i C m-)]- After some algebra, this 
ends up being equivalent to the condition £ < £ T (t). Physically, this condition simply arises 
from the time required for a loop of size £ to shrink down to the vorton length £ v . The 
quantity £ T (t) is also the upper bound of the vorton length spectrum at any time, their size 
ranging from A to £ T (t). 

Concerning the first line of Eq. (3.5), accounting for proto- vortons and doomed loops, 
the integral over can again be performed explicitly owing to the Dirac function. However, 
one has to distinguish two cases according to the value of Nq = [£/X + 7d(i — tcur) /A] 1 / 2 , 
the zero of the Dirac function argument. Either Nq < and £ v (Nq) = 0, or Nq > A* and 
^v(Ao) = Nq/VU. Again, these conditions can be recast in terms of lengths by defining a 
new length scale 

= A [A 2 - 7d (t - t cur )] , (3.10) 

which is a decreasing function of time. At any time, all loops having £ < £*(t) are doomed, 
they are the ones associated with N < A* and will disappear by decay. For those, we have 
set above £ v = and the Heaviside function in the first line of Eq. (3.5) equals unity. All the 
others loops, having £ > £*(t) are proto- vortons, i.e. in the shrinking stage before becoming 
vortons. This can be explicitly seen by rewriting the Heaviside function as 



&[£ - £ V (N )] = G[4(i) -£] + &[£ - £*(t)} @[£ - £ T (t)] . 
Combining all terms together, the loop distribution function finally reads 

a 3 (t) — (t, £) = af ni {9[4(t) -£] + &[£ - *•(*)] &[£ - e T (t)]} + 7d(< - km)} 



(3.11) 



2al i N^e[£ T (t)-£]Q(£-X)M 1 



N 



A 



+ 7d(^c 



ti 



(3.12) 



where A* is given by Eq. (3.4), £ T (t) and £*(t) being defined in Eqs. (3.9) and (3.10), respec- 
tively. As we have just discussed, the second line of this expression accounts for all vortons 
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present at time t while the first describe both proto-vortons and doomed loops. Let us stress 
that the doomed loops, i.e. the domain £ < £*(t) exists only during a transient period after 
which they completely disappear. This happens at the time t = t* solution of = 0, i.e. 

for 

t* = tcur H ■ (3.13) 

7d 

The loop distribution of Eq. (3.12) generalizes the approach of Ref. [64] for any initial loop dis- 
tributions while tracking at all times both the proto-vorton and vorton populations. Taking 
an initial loop distribution peaked at a particular length gives back the results of Refs. [64, 65] 
in the asymptotic limit t S> i* and if we neglect the proto-vortons. As a result, we do not 
expect significant new constraints to be derived from Eq. (3.12), i.e. in the case of a vanishing 
loop production function. However, as can be seen from Eq. (2.43), for V(£,t) ^ 0, as it is 
found in the numerical simulations, the vorton population should be significantly enhanced. 
We let however for a future work the derivation of the associated cosmological constraints. 

4 Conclusion 

We have presented a new way of treating the evolution of a string network, taking into 
account in principle any kind of interaction between loops and long strings, and in particular 
the appearance of an ensemble of vortons or quasi- vortons [64], i.e. loop configurations whose 
internal structure induces a different, possibly vanishing, decay rate. These vortons could 
plague GUT models having a cosmic string network building as a consequence of primordial 
symmetry breaking, i.e., basically all cosmologically-compatible GUT [74, 75]. Indeed, in 
realistic models, the strings that form are expected to couple to various fields [51], some 
of which leading to current-carrying vortices [53] in such a way that the strings are no 
longer Nambu-Goto like as usually assumed [67, 76-78]. This leads to the formation of very- 
slowly decaying, or even stable, matter scaling objects dubbed vortons, that can change the 
subsequent evolution in a drastic way. Up to now, only rough evaluation of their contribution 
has been proposed using their initial distribution function. 

In this work we have gone one step further into understanding in more details the vorton 
distribution and the relevant cosmological constraints. We have presented a new Boltzmann 
approach governing the evolution of the number density n(t, £, N) of current carrying cosmic 
string loops. Under some assumptions, we have been able to find an explicit solution starting 
from any initial distributions and for any given loop production functions: the most general 
results are given by Eqs. (2.24), (2.37) and (2.39). 

We have also shown how this method was extending previous results on the cosmologi- 
cal evolution of vortons in the absence of loop production [64, 65], while remaining perfectly 
compatible in the appropriate limits. Our results could be readily applied to the loop pro- 
duction function found in Ref. [31, 35, 36] to revisit the vorton constraints in presence of 
loop production, as should be done in a forthcoming work. Another possible extension could 
be the determination of a complete analytic loop production function but from a system of 
Boltzmann equations that would explicitly couple long strings and loops with collision and 
fragmentation terms. Indeed, much still deserves to be done in this area: even though cosmic 
strings have long ago been ruled out as the main source of density perturbations, it does not 
mean that they ought not be used as a very powerful tool to provide constraints on theories 
otherwise unreachable. 
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